EMPG-03-10 



Auxiliary matrices for the six-vertex model at q = 1 II. 
Bethe roots, complete strings and the Drinfeld polynomial 



Christian KorfT 

School of Mathematics, University of Edinburgh 
May field Road, Edinburgh EH9 3JZ, U.K. 



Abstract 

The spectra of recently constructed auxiliary matrices for the six-vertex model 
respectively the spin s = 1/2 Heisenberg chain at roots of unity are investigated. 
Two conjectures are formulated both of which are proven for N = 3 and are ver- 
ified numerically for several examples with N > 3. The first conjecture identifies 
an abelian subset of auxiliary matrices whose eigenvalues are polynomials in the 
spectral variable. The zeroes of these polynomials are shown to fall into two sets. 
One consists of the solutions to the Bethe ansatz equations which determine the 
eigenvalues of the six-vertex transfer matrix. The other set of zeroes contains the 
complete strings which encode the information on the degeneracies of the model due 
to the loop symmetry s?2 present at roots of 1. The second conjecture then states 
a polynomial identity which relates the complete string centres to the Bethe roots 
allowing one to determine the dimension of the degenerate eigenspaces. Its proof for 
N = 3 involves the derivation of a new functional equation for the auxiliary matrices 
and the six-vertex transfer matrix. Moreover, it is demonstrated in several explicit 
examples that the complete strings coincide with the classical analogue of the Drin- 
feld polynomial. The latter is used to classify the finite-dimensional irreducible 
representations of the loop algebra s^- This suggests that the constructed auxil- 
iary matrices not only enable one to solve the six-vertex model but also completely 
characterize the decomposition of its eigenspaces into irreducible representations of 
the underlying loop symmetry. 
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1 Introduction 



This article is a continuation of a previous work on auxiliary matrices for the six- 
vertex model respectively the s = 1/2 Heisenberg chain at roots of 1, 

M _! 

17 \ ^ _i , _« v + ( z „z i\ x,y,z _ x,y,z \ 

H = 2^ °m°m+l + ff m a m+\ H ^ K^m^m+l ~ *■)■> °M+1 = 1 ■ I 1 ) 

m=l 

Here <Jm V ' z denote the respective Pauli matrices acting on the m th site. Throughout 
this article the deformation parameter q will be assumed to be a primitive N th root of 
unity with N > 3. While the first part of this work was mainly concerned with the 
construction of the auxiliary matrices and their geometric structure, the present paper 
focuses on the nature of their spectra. This will allow us to solve the eigenvalue problem 
of the six- vertex transfer matrix and thus the Hamiltonian <[T]>. In particular, we will 
derive the Bethe ansatz equations [2 IU 03 

sinhi(nf +i 1 ) \ M = |* sinh^f -nf + 2i 7 ) ^ = ^ 
sinh \{u® — i'y) J sinh \{u? — uf — 2ij) ' 

from representation theory and identify their solutions as zeroes of the auxiliary matrices' 
eigenvalues. Furthermore, we will demonstrate in concrete examples how the spectra of 
the auxiliary matrices yield information about the degeneracies connected with the loop 
symmetry sl% present at roots of 1 [HUZj. Before we describe the results of this paper 
in detail we give a short introduction into the method of auxiliary matrices and a brief 
overview of the results obtained in pQ. 

1.1 Auxiliary matrices from quantum group theory 

The concept of auxiliary matrices was originally introduced by Baxter in the context 
of his solution to the eight-vertex model [HI EEB CCD EE2| and motivated by the lack 
of spin-conservation. His approach is described in detail in ^Hj. While the six- vertex 
model preserves the total spin, its infinite-dimensional symmetry algebra sfa at roots 
of unity does not. It is for this reason that auxiliary matrices provide the appropriate 
approach for the discussion of the spectrum and the degenerate eigenspaces. Unlike 
the coordinate space Bethe ansatz Jl] they do not rely on spin-conservation. Also 
the algebraic Bethe ansatz ^1 nas serious deficiencies. Away from a root of unity 
the entries of the monodromy matrix provide a spectrum generating algebra providing a 
complete set of eigenstates. This ceases to be true when q N = 1 as then certain operator 
products in the Yang-Baxter algebra vanish. If one wants to resolve the structure of the 
degenerate eigenspaces at roots of unity the concept of auxiliary matrices is therefore 
the only method left. 

However, as far as the six-vertex model is concerned an explicit construction for an 
auxiliary matrix has only been given by Baxter for the sectors of vanishing total spin, 
cf. formula (101) in [23- His expression applies to generic values of the deformation 
parameter. More recently, auxiliary matrices related to the six- vertex model for q N ^ 1 
have been studied using infinite-dimensional representations of quantum groups ^01 EEE1 
I17U19| . In |19j the connection with the lattice model has been explicitly investigated and 
Baxter's result has been extended to spin sectors different from zero involving formal 
infinite power series. 
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Quantum group theory also played the key role in the construction of auxiliary 
matrices at q N = 1 discussed in [Q. It needs to be emphasized that in contrast to 
[W\ \1H\ I17| I19j the construction in pTJ uses finite-dimensional representations which are 
special to the root-of-unity case. Furthermore, this approach differs in several essential 
points from the one outlined by Baxter, see Sections 1.2 and 1.3 in [Q. Nevertheless, 
the key idea remains the same. 

Recall that the statistical lattice model is defined in terms of the transfer matrix T(z) 
(defined in equation (|2()[) below) which besides the deformation parameter q depends on a 
spectral variable z. The spin-chain Hamiltonian (^TJ) is obtained from the transfer matrix 
by taking the logarithmic derivative with respect to z and setting z = 1 afterwards. One 
now introduces an additional matrix Q which following Baxter is called "auxiliary" . Its 
defining property is the solution of a suitable functional relation which allows one to 
solve the eigenvalue problem of the transfer matrix T in terms of Q. The main result of 

was to explicitly solve the following operator functional equation, 

Q p (z)T(z) = Mz) M Q P >(zq 2 ) + Mz) M Qp"(zq~ 2 ) ■ (3) 

Here 4>i , 4> 2 are scalar functions (cf. equation (|37|) in this article) . The auxiliary matrices 
Qp-, Qp'j Qp" depend on additional complex parameters p = (x, y, z, c = fi + /i _1 ) G C 4 
with z / 1, whose appearance is connected with the enhanced symmetry of the six- vertex 
model at roots of 1. These parameters define points on the following three-dimensional 
complex hypersurface [20l I2TJ l22 | 123] 

ci rv n' i -i\ n' , -TV' at/ f N, if N odd , . 

SpecZ: xy^(z + z 1 )=/ + ^ ) N = { N/2, if N even ' ^ 

The points p' = (x, q N 'y, q N z, fiq + p>~ 1 q~ 1 ) and p" = (x, q N y, q N 'z, pq~ x + H~ 1 q) in 
the functional equation (j3J) are determined by representation theory. When vV is even 
one has to make the further restriction x = y = 0. For details and the explicit definition 
of the auxiliary matrices we refer the reader to T_. For a subvariety of SpecZ their 
definition is given below (see equation IJ2fijl ) in order to keep this article self-contained. 
All matrices in the functional equation © have been proven to commute with each 
other whence it can be written in terms of eigenvalues. 

There are two main advantages of considering the auxiliary matrix Q p as the central 
object instead of the transfer matrix. First, the Bethe roots solving (j2j) can be 
directly obtained as zeroes of the auxiliary matrices' eigenvalues 

= Q p (zf)T(zf) = Mzf) M Q P izfq 2 ) + <t> 2 (zf) M Q P »(zf<l- 2 ), 4 = <?* ' q' 1 • (5) 

Second and more importantly, Q p breaks the infinite-dimensional symmetry of the six- 
vertex model at roots of unity and therefore is in general non-degenerate. Employing 
the functional equation (JHJ) one can show that this implies that the eigenvalues of the 
auxiliary matrices must contain factors which are (^-periodic. Consequently, the eigen- 
values of Q p possess additional zeroes besides the Bethe roots which are called complete 
strings ^21 12] 

q N = l: Q p (zfq 2e )=0, £ = 0, 1, 2, .... N' - 1 . (6) 

The occurrence of these complete strings at finite length M of the spin-chain is charac- 
teristic to the root-of-unity case. Note that in contrast to the Bethe roots z? the string 
centre z s is not determined via the functional equation © . This freedom is at the heart 
of understanding the infinite-dimensional symmetry at roots of unity. As we will see in 
this article the structure of the degenerate eigenspaces of the transfer matrix and the 
Hamiltonian Q is completely described by the complete strings. 
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1.2 Results and outline of the article 



The investigation of the spectra of the auxiliary matrices Q p for arbitrary points on the 
hypersurface Q is a quite complicated task since they form a non-abelian set. That 
is, for a generic pair p\,p2 € SpecZ and any pair of spectral variables z,w £ C the 
corresponding auxiliary matrices do in general not commute. Instead one finds that 
in order to ensure [Q pi (z), Q P2 (w)] = one has to enforce the relations 

—eN' —eN' eN' eN' 

xiz £JV x 2 w yi-z y2^ - - o l (7) 



1 — Zi 1 - z 2 ' 1 - z, 1 1 



-1 ' 



Note that this is sufficient to guarantee that all operators in Q commute. While the 
non-abelian character of the auxiliary matrices makes them more powerful as a symmetry 
it complicates the calculation of the eigenvalues, since the eigenvectors may depend on 
the additional parameters p, the spectral variable z and the deformation parameter q. 
We will therefore focus only on an abelian subset for which the eigenvectors exclusively 
depend on q and the spectra consist of polynomials in the spectral variable z. This 
abelian subset is defined in the following conjecture which we will prove for the case 
N = 3 in this article. 

CONJECTURE 1. For integer N > 3 consider the subvariety in the hypersurface ^ 
defined by 

Vil = (0,0, fi~ N ' , fi + Li' 1 ) G Spec Z, /i£C x . (8) 

Denote the corresponding auxiliary matrices by Qfi(z) = Q P A Z )- An explicit definition 
will be given below, cf. equations \2b]) . 128\) and \29\) . Then one has the commutation 
relations 

[Q^z),Q v (w))=0, (i,!/eC x , z,!ueC. (9) 

Because of this relation we refer to the one-parameter subset of auxiliary matrices 
as "abelian". 

For N = 3 this conjecture will be shown to be valid by explicitly constructing 
the intertwiner of the quantum loop algebra U q {sl2) associated with the evaluation 
representations 71^,71^" (see definitions (|29j). (|E3*|) in the text and (|8TI) in the appendix). 
That is, we will demonstrate for q 3 = 1 that the tensor products t^z 1 " and 7rSr (8>7t^ m 

are isomorphic. Numerical calculations have also been performed for N = 4, 5, 6, 8 and 
the conjecture has been found to be valid. 

The conjecture that this assertion holds true for all N is motivated by the observation 
that the necessary criteria for the existence of such an intertwiner - which are identical 
with the one shown in - are satisfied. These necessary criteria turn out to be 
sufficient for the existence when Xj, 7^ 0. The corresponding intertwiners have first 
been obtained in When Xj = yj = the above criteria are obviously trivially 

satisfied for any values of z,w and zi,Z2- Unfortunately, the parametrization used in 
|25| does not allow one to take the limit Xj,yj — ► and to obtain the corresponding 
intertwiners for the subvariety Thus, we need to prove the existence for this special 
case which is done by explicit construction in the appendix of this article for N = 3. 
The case of arbitrary N is left to future work. 

It follows from their definition given in equation l|26|) below that each auxiliary matrix 
can be decomposed as 

M 
m=0 
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where the coefficients are independent of the spectral variable z. If Q holds true 
all the coefficients commute. In fact, one has [Qj^\Q^] = 0. Hence, every eigenvalue 
of the auxiliary matrix Qu(z) can be written in terms of polynomials with the most 
general form being 

Q„{z) = N^z n °°P B (z)P^z)P s (z N \u) • (11) 

Here A/" M = N^(q) is a normalization constant not depending on z and the three poly- 
nomials are given by 

TLB 

Pb(z) = \[{z-zf), (12) 

i=i 

n 

P*{*) = Uiz-zjifi)), (13) 
i=i 

ns ns 

p S (z N '^) = n n (*- = - ■ 

j=l 1& N , j=l 

The first polynomial Pb contains the zeroes zf = z?(q) of the eigenvalue which do not 
depend on the parameter \x and for which there is at least one t G Z^v' such that z^q 2e is 
not a zero of Pg. We will show below that these zeroes are finite solutions of the Bethe 
ansatz equations ([%]). whence the notation. Via the polynomial Pb(z) they determine 
(up to a possible sign factor) the eigenvalue of the six-vertex transfer matrix associated 
with JTTJ), 

T(z) = M*) M * n ~?0± + M*) M <l- n °° • (15) 

The power of the monomial in 1)11)1 gives the number of the "Bethe roots at infinity" , 
see e.g. [21] and references therein. This expression refers to the parametrization in © 
with zf = e u > q- 1 -> 0. The appearance of "infinite" Bethe roots is another feature 
characteristic to the model at roots of 1. It signals the breakdown of the familiar 
formula that the number ns of Bethe roots is related to the number of down spins in 
the corresponding eigenstate. 

The second polynomial P^ accounts for the possibility of zeroes depending on q and 
fi. They occur because the additional parameters p shift in the functional equation 
(j2J). Again we allow only for zeroes Zj(fi) for which there is at least one £ E Z^r/ such 
that Zj(fi)q 2e is not a zero of P„. Using the transformation behaviour of the auxiliary 
matrices under spin-reversal we will show that 

n B 

Pa(z) = H(Z - zffl 2 ) = ^P B (Z^ 2 ) . (16) 
J=l 

The last factor P$ contains the contribution of the complete iV'-strings where 
the string centre Zj may or may not depend on fj,. The additional dependence on the 
deformation parameter q is suppressed in the notation. As mentioned before the contri- 
bution of the complete strings encodes the information on the degeneracies connected 
with the loop symmetry at roots of unity. In particular, the number ns of complete 
strings is related to the dimension of the corresponding degenerate eigenspace of the 
transfer matrix. In order to arrive at this result we need to determine the number of 
possible eigenvalues (jTTj) in a degenerate eigenspace of the transfer matrix. 
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In the case N = 3 this information will be derived from the following functional 
equation, 

N = 3 : Q^z)Q u (zv 2 q 2 ) = Q^ q (zv 2 q 2 ) [q M (z - l) M T(zq) + (zq 2 - 1) M ] . (17) 

The above identity relates a product of two complete string contributions (|14j) to a single 
one, thus imposing severe restrictions on the possible /U-dependence of the string centres 
Zj(fi). One finds that only two possibilities are allowed: namely, one has 

either zf(ji) = zf or zf(ji) = z^ 2 . (18) 

Here Zj denotes the (constant) value of the string centres in the limit fj, — ► 1. These 
values are fixed in terms of the Bethe roots via the following remarkable identity which 
is also deduced from (|17[). 

. n 2(£+l) nao ( 2i_-,\M 

, „) - £_ ■ (10) 

Both results taken together imply that for each eigenvalue Q15j) of the transfer matrix 
which possesses a fixed number of Bethe roots, finite and infinite ones, there are 2 ns 
possible eigenvalues of the auxiliary matrix. Since the auxiliary matrices break the 
infinite-dimensional symmetry of the six-vertex model as well as spin-reversal symmetry 
they are non-degenerate. Thus, the corresponding 2 ns eigenstates yield a basis for the 
degenerate eigenspace of the transfer matrix. 

Note that the crucial functional equation (j!7j) is only valid for N = 3 and must be 
modified for N > 3. Nevertheless, the outcome ought to hold true in general leading to 
the formulation of the second conjecture. 

CONJECTURE 2. The identity hiy\) and the restriction U8\) on the [i- dependence 
of the complete string centres not only holds true for N = 3 but applies to all primitive 
roots of unity of order N > 3. This in particular implies that each eigenvalue of the 
six-vertex transfer matrix allows for 2 ns eigenvalues of the auxiliary matrix. 

While we will provide a proof of this identity only for iV = 3we performed numerical 
checks for N = 5, 6, 8 verifying the second conjecture also in these cases. Note also that 
this assertion coincides with previous results in the literature [HJ 1241 1281 127j obtained 
by numercial calculations and use of the loop algebra symmetry s?2 which has been 
established in the commensurate sectors 2S Z = OmodiV jSl[7j. In this context Fabricius 
and McCoy suggested in [2Sj an expression for the classical analogue of the Drinfeld 
polynomial [2JJ|. The latter describes the finite-dimensional irreducible representations 
of the loop algebra [30] . 

The main result of this article is that the contribution of the complete strings (|19|) 
of the auxiliary matrices constructed in [Q coincides with the proposed expression for 
the classical analogue of the Drinfeld polynomial. We will address this point in the last 
section of this paper when discussing concrete examples. They suggest that the spectra 
of the auxiliary matrices describe the decomposition of the eigenspaces into irreducible 
representations of the loop algebra. This is of particular importance as the auxiliary 
matrices have been defined for all spin sectors, while the loop algebra generators have 
only been constructed for the sectors where the total spin is a multiple of the order N. 

The outline of the article is as follows. In Section 2 we introduce our conventions in 
defining the six-vertex model and review the definition and properties of the auxiliary 
matrices. Section 3 discusses the implications of the functional equation (j3J) for the form 
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of the eigenvalues Section 4 deals with the transformation under spin-reversal, a 

symmetry which is broken by the auxiliary matrices. Section 5 shows that the eigenvalues 
of the auxiliary matrices occur always in pairs of opposite momenta and gives the relation 
between them. Section 6 contains the proof of the crucial functional equation 1)17(1 . 
Section 7 summarizes the results. In particular, the connection with the representation 
theory of the loop algebra s?2 is made. 

2 Definitions 

In order to keep this paper self-contained we briefly recall our conventions for the defi- 
nition of the six-vertex model. The transfer matrix is given as the following trace over 
an operator product 

T(z) = trRo M (z)R M-i(z) ■ ■ ■ R i(z) (20) 

with 

= — ^ — 1 <2> 1 H — a z ® a z + ca + <g> a + c a ®a + , a ± = (21) 

being defined over C 2 (g> C 2 in terms of the Boltzmann weights of the six-allowed vertex 
configurations, 

(l-z)q 1-q 2 , 

a = p, b = p- T , c = p- n ' c=cz. (22) 

1 — zq A 1 — zq A 

For convenience we shall henceforth set the arbitrary normalization factor to p = 1 . The 
lower indices in (|20|) indicate on which pair of spaces the iZ-matrix acts in the (M + 1)- 
fold tensor product of C 2 . The explicit dependence on the parameter q is suppressed in 
the notation. 

The six-vertex transfer matrix possesses a number of finite symmetries given by the 
vanishing of the following commutators 

[T(z),S z ] = [T(z),m] = [T(z),e\ = , (23) 

where the respective operators are defined as follows 

1 M 

S z = -^a z m , m = a x ®---®a x , 6 = a z (8) • • • (8) a z = (_i)^/2-|S*l . ( 2 4) 

m=l 

The first operator is the total spin, the second invokes spin-reversal and the third has 
eigenvalue +1 or —1 depending whether the number of down spins n in a state is even 
or odd. The finite symmetries and the properties of the Boltzmann weights (|22|) can be 
used to derive for spin-chains of even length, M = 2M' , the following useful relations of 
the transfer matrix, 

T{z,q~ l ) = T(z-\q), 
T(z,-q) = eT(z,q) = T(z,q)e, 
T(zq-\q) = biz-^Tiz-^qf . (25) 

Here we have temporarily introduced the explicit dependence on the deformation pa- 
rameter q in the notation. The (|25|) transformation properties impose restrictions on 
the spectrum of the transfer matrix. Henceforth, we shall assume M to be even. 
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2.1 A one parameter family of auxiliary matrices 

As explained in the introduction only a subclass of the auxiliary matrices constructed 
in will be considered, namely the ones associated with nilpotent representations. In 
terms of the hypersurface Q this class of auxiliary matrices corresponds to the points 
jSJ). We now explicitly define this one-parameter family of auxiliary matrices setting 

Quiz) = \xL* M {z/n)I% M _ x {z/n) ■ ■ ■ L^z/p,), L» m G End(Vo ® V m ), /i G C x . (26) 

Here the operators in the trace can be expressed as 2 x 2 matrices with operator-valued 
entries 

U t =(^j t = A^®a + <j- + B^®a + + C^®<j- +D^®<j-ct + . (27) 

The operators A^, 13^,0^, G End(C Ar ') are given as 

A„(w) = wq^(t)-^{t)-\ 

C, = (g-g-yieKli)- 1 , 
D^w) = wq^(ty l -^{t) (28) 

with the N' x N' matrices ^(t), 7r^(e), Tr^(f) defined through the following action on 
the standard basis {v n } in C^' (13 EU * , 



^{tfv n = fl^q-^Vn, 7T^f) Vn = V n+l , 1^U>N'-1 = 

H + H-^-luF-H-iq- 2 " 



(q 



-1\2 



V n -1 ■ (29) 



The matrices ^(i) 2 , 7T^(e), 7r^(/) define a representation 7r^ of the quantum group 
U q {sl2) at = 1, i.e. they are subject to the relations 

Tr^TT^e)^)- 1 = gTr^e), 
7r^(i)^(/)^(i)- 1 = g-V(t), 

[^),^(/)] = ^ ft) ;_-;_? r2 ■ (30) 
The elements generating the centre of the quantum group take the values 

n ^ff = ^(ef = 0, ^) 27V ' =^ N ', (31) 

and 

+ g -l^( t )-2 + (g - (? - 1 ) 2 ^(/)^ (e ) = ^ + M "l . (32) 

There are several advantages of defining the auxiliary matrix in terms of the repre- 
sentation ir^. As a consequence of the quantum group relations one has the identity 

L» 12 (w/z)L^(w)R 23 (z) = R 23 (z)L^(w)L^( W /z) (33) 



"This parametrization of the auxiliary matrices and the root-of- unity representation slightly differs 
from the one used in |2()II21) respectively [T|. Instead of using the parameter A (cf equation (43), Section 
2 m CP) h is more favourable to use the variable /i = \~ 1 q~ 1 as this facilitates the identification when 
taking the nilpotent limit from a generic cyclic representation p — <p(£,C, A). 
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which implies that the auxiliary matrix and the transfer matrix commute, 

[Q lx (w),T(z)]=0. (34) 

In addition, one derives from the following non-split exact sequence of evaluation rep- 
resentations 7Tw 

0->7r$^-><®7r°-> TT^" 1 -► 0, w = w'q" 1 = w"q = z/fi (35) 
the specialization of the functional equation (j3J) to the subvariety |JHJ) 

Q^(z)T(z) = Mz) M Q M (zq 2 ) + ^(^) M Q W -.(C 2 ) • (36) 

Here ir® is the root of unity limit of the two-dimensional fundamental representation 
defining the six- vertex model. The corresponding point on Spec Z is p° = (0, 0, q N ,q 2 + 
q~ 2 ). The complex numbers w,z play the role of evaluation parameters respectively 
spectral variables. The coefficient functions are given by* 

<t>i{z,q) = b(z,q)q~% and <f> 2 {z) = q* . (37) 

All matrices in the equation (|36j) have been shown to commute, whence they can be 
simultaneously diagonalised and the eigenvalues of T can be expressed in terms of the 
eigenvalues of the respective auxiliary matrices. 

For JV = 3 we will show that the auxiliary matrices obey the stronger commutation 
relations ©• For N > 3 we assume that Conjecture 1 holds for the reasons stated in 
the introduction. 

Recall from that the auxiliary matrices for generic p £ Spec Z break all the finite 
symmetries of the six-vertex model. However, the one-parameter family (|26|) which 
constitutes a subvariety preserves two of the finite symmetries of the six- vertex transfer 
matrix, namely 

[Q ll (z),S z ] = [Q ll (z),&}=0 . (38) 
Spin-reversal symmetry on the other hand remains broken 11, 

K Q^z)m = Qp-i (z^ 2 ) = (-zq/fi) M Q^q-Yf ■ (39) 

Employing the U q (sl2) algebra automorphism e — ► /, / — ► e, t ±l — ► t ±l , q — > q^ 1 one 
proves the additional relation 

Q„(z,q)=Q f ,(zq 2 ,q- 1 ) t (40) 
which allows one to derive the adjoint of the auxiliary matrix 

Q?(z, Q)* = Quiz, I' 1 ) 1 = Qn(zq~ 2 ,q) ■ (41) 

These identities will prove crucial in the following investigation of the eigenvalues (jlljl - 

^These coefficients are obtained from the one in p@ by setting p±(z) = (zq)^~ in equation (69), 
Section 3. 
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3 The TQ equation in terms of eigenvalues 



We start our analysis of the spectra of the auxiliary matrices by inserting the expression 
(jlljl in the functional equation (|3(ij) . Recall from the introduction that (jlljl is the most 
general form of the eigenvalues provided © is true. By abuse of notation we will denote 
the operators and their eigenvalues by the same symbols. We obtain from IfMfijl . 

T(z)N^P B (z)P fl (z)P s (z NI \n) = 

Uz) M M m q 2n -P B {zq 2 )P, q (zq 2 )P s {z N '^q) 

+ faW^M-i PB(zq- 2 )P M -i(zq- 2 )P S (z N ' ^q- 1 ) 

The eigenvalues of the transfer matrix Q2U|) are independent of /i, whence the polynomials 
Pfii Pfiq, P^q- 1 must cancel on both sides of the functional equation (|36|) except for a 
constant factor q ±2n . Up to a possible renumeration of the zeroes Zj(/i), this implies the 
following relation 

Zj(nq) =Zj{fi)q 2 . (42) 

By the same argument one deduces that the complete string centres Zj(fx) can only 
differ by powers of q for \i — > /xq ±:L , 

Ps(z N ',fi)=P s (z N ', m ), => P s (z N ',») = P s (z N ',ti N ). (43) 

In addition, we can conclude that the ratios A/^/A/^, A^-i/A^ of the normalization 
factors are independent of fi, 

M m /N» = A^/AT M -i . (44) 
Thus, we deduce the following preliminary form of the eigenvalues of the transfer matrix 



T{z) = Mz) M Q 2n ~ +2nB ^^j^ + Uz?U~ 2n - 2nB ^ ^f^f ■ (45) 

This form is "preliminary" as we will determine the ratios of the normalization constants 
below. Evaluating the left-hand-side of the TQ-relation at a zero z = z? we obtain also 
a preliminary form of the Bethe ansatz equations, 

= M*?) M Mn q 2n °°P B {zfq 2 ) + <j> 2 {zf) M q- 2n °°N m -i P B {zfq- 2 ) . (46) 

These equations ensure that the eigenvalues of the transfer matrix have residue zero 
when the limit z — > z? is taken. Below we will see that the zeroes of the polynomial 
Pb in fact coincide with the finite solutions of the Bethe ansatz equations ©, i.e. they 
are the Bethe roots at q N = 1. Note that it might also happen that the zeroes z? = 1 
and z? = q~ 2 simultaneously occur. While this possibility looks problematic in light of 
the parametrization used in (j2J, equation (|4Hj) shows that it does not pose a problem as 
both sides of the equation then vanish. Also the limit z — * 1 yielding the momentum 
eigenvalue in (|45|) stays well-defined. We therefore include these zeroes in the set of 
Bethe roots. 
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4 Transformation under spin-reversal 

We now discuss the implications of spin-reversal. First note that we can deduce from 
(jSU, © an d the integrability of the six- vertex model, [T(z),T(w)] = 0, that the auxil- 
iary and transfer matrix have common eigenvectors which neither depend on the spectral 
variable z nor on the parameter [i. Let v = v(q) be such a common eigenvector with 
eigenvalue (jllj) . Then we find according to Q39JI that 

-2\ 



jri 00 .,-2n 00 r> /„„-2\n /~,.-2\ D f„N',.-2N' ,,-1 



A^-i z^ti-^PBiz^P^iz^Psiz" vT l )9tv . (47) 



Furthermore, the above eigenvalue must satisfy the functional relation (j36j) . 

Q^{z)T{z)y\v = [Mz) M Q m (zq 2 ) + Mz) M Q m -i(zq~ 2 )]Kv 

VKQ^(zn~ 2 )T{z)v = ^[0 1 (z) M Q At - l9 - 1 (z/i- 2 )+0 2 (z) M Q^ 1 ,(z^- 2 )] ? ;. (48) 

Here we have exploited that the transfer matrix is invariant under spin-reversal. In the 
previous section we verified that the eigenvalues of T[z) have 3ns zeroes at z^,z^q ±2 . 
These zeroes correspond to the finite Bethe roots as we will see shortly. In equation (|48|) 
these zeroes must originate from the polynomial 

n 
3=1 

as the factors Pb(z[i~ 2 ) and the contribution of the complete strings Ps(z N ' fi~~ 2N ' , fi~ N ) 
cancel on both sides of the equality sign. In fact, replacing fi — > we obtain the 
expression 

We saw already earlier that T(z) does not contain any poles, i.e. the zeroes Zj([i)fi~ 2 
have now to satisfy the preliminary Bethe ansatz equations, 

Consequently, we have the zeroes Zj{^i)ix~ 2 , Zj(fj,)/j,~ 2 q ±2 of the eigenvalue which must 
coincide with the zeroes z? , z^q^ 2 . Hence, after a possible renumeration we are lead to 
the conclusion 

zf = z j (»)»- 2 (49) 

which proves relation (|16|) . Using this identification we can now determine the ratios 
of the normalization constants by comparing with the earlier expression (|45|) for the 
transfer matrix eigenvalue and employing (|44|). 



fc) 2 = = ^+4n B ^ ^ = ±Kq ~n^2n B (5Q) 



Inserting this result back into (|45|) the final expression (|15j) for the six-vertex transfer 
matrix eigenvalue associated with (|11|) is obtained up to a possible sign factor. This 
ambiguity is due to the square root in (|50|). 
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The missing step in order to compare this result with the outcome of the coordinate 
space Bethe ansatz is to verify that the equations (|46[) coincide with ©. Employing 
ijllj) a simple calculation gives 



V .B„2\ M n B R„2 



Z 3 Q \ _ g 2n 00 +2n B -M TT 2 / ( l 



Except for the additional phase in front of the product this coincides with (J2J) if we 
identify 2 = e u q , q = e* 7 . For real eigenvectors we will show momentarily that the 
phase factor is equal to one. In the general case of complex eigenvectors we do not 
have a proof that the phase factor is always trivial. However, numerical calculations for 
N = 3,4,5,6 and spin chains up to the length M = 11 have so far not produced any 
counter example. 



5 Pairs of eigenvalues 

In this section we are going to exploit the relations (|39j) and (|4()[) of the auxiliary matrix 
to show that the eigenvalues occur in pairs. Combining these two identities we obtain 

Q^Z,q) = (-zq/^Q^iz-'q-^q) 1 = (-zq/ » M Q^{z~\q- 1 ) . (52) 

According to @ we can find eigenvectors v which only depend on the deformation 
parameter, i.e. v = v(q). Equation (|52jl then shows that the eigenvectors come either 
in pairs, (v(q), v(q~ 1 )) with 

Q fl (z,q)v(q)=Ar fl (q)z n -P B (z,q)P^z,q)P s (z N ',fX N )v(q) , (53) 

and 

Qix{z,q)v(q~ x ) = 

(-zq/fl) M K-i (q- 1 ) z- n °°P B (z-\ q- l )P„-r {z~\ q~ ^Ps^"' , M^M?" 1 ), 



or are real, v{q) = f(g _1 ) = v(q), in case of which the two eigenvalues above must be 
equal. From (|25|) we infer that the two eigenvectors (v(q), v(q~ 1 )) have momentum k of 
opposite sign with 



Ak 



^ff 1 -^ 



Rewriting the eigenvalue of v(q 1 ) in the form one deduces that the number of 
infinite Bethe roots transforms according to 



M- noo -2n B - n s N', (54) 



while the finite Bethe roots and complete string centres of the respective eigenvalues are 
related by the transformations 

zf (g) -+ 1/zf (g- 1 ) and zf (ji N ) -+ 1/zf (^) . (55) 

Finally one finds for the normalization constants the mapping 

n B n s 
3=1 3=1 
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The case of real eigenvectors can only happen when the corresponding eigenvalues 
of the translation operator T(z = l,q) are real, i.e. k = 0,ir. Then 1)52(1 becomes 
an identity in terms of eigenvalues and Bethe roots as well as string centres must be 
invariant under the transformation laws (|54j) . (|55|) . This implies for v (q) = v(q~ 1 ) the 
identities 

M = 2noo + 2n B + n s N', (57) 

IB 

P B {z-\q- l ) = {-z)-^P B (z,q)J{zf{q)-\ (58) 

3=1 

ns 

P s (z- N ',»- N ) = P s (z N ',» N )(-z N ')-^llzf (59) 

3=1 

Note that (|57|) fixes the phase factor in l|51|) to be q~ n s N ' = (— i) n s(N+i) anc [ completes 
the derivation of the Bethe ansatz equations ((2j) for odd roots of unity and real eigen- 
vectors. For N even we will see momentarily that there is always an even number of 
complete strings, whence (j2j) also applies in this case. 

The relation (|57|) also implies that in a degenerate eigenspace of the transfer matrix 
with real eigenvectors and a fixed set of Bethe roots the number of complete strings is 
constant. Below we will see for N = 3 that this also holds true for complex eigenvectors 
by proving the functional equation (jTT|) and the identity (fT§|) . 

Inserting the expressions (|58|) and (|59|) into the identity (|52|) yields the following 
equation for the complete string centres, 

wM^ M+2n ° n z 3^) N> = f nB im) 2 = 1 • (6°) 

3=1 3=1 

Here we have used that 

riB 

limT(M) = (-q)- nB U z f ^ = ±X ■ ( 61 ) 

3=1 

Exploiting that Zj (p~ 1 q ±l ) N ' = ^(/i -1 )^' and (|5U() leads to the further restriction 

( M ) M+2ns (AA M /AA M -i g -i) = g M - 2nfl - 2n »M Jlf+2nfl (^/^ t - 1 ) ^ ^ = 1 • ( 62 ) 

Thus, in the case of even roots of unity there is always an even number of complete 
strings. This completes also the derivation of the Bethe ansatz equations ((21) for even 
roots of unity by showing that the phase factor in (|5Tj) is equal to one. 

Note that up to this point all relations have been derived for general N' > 3. Thus, 
while the conjecture © is only proven for N = 3 in this article, the derivation of the 
spectrum of the auxiliary matrices applies to all roots of unity once the commutation 
relations © are established. 

6 A new functional equation for N = 3 

In this section we prove for g 3 = 1 the two important formulas (|19|) and (|18|) employing 
the functional equation (|17|h While the explicit form of this functional equation is 
characteristic to the case N = 3 the representation theoretic method applied to derive it 
is not. Indeed, the line of argument which employs the decomposition of tensor products 
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of evaluation representations via exact sequences also applies to the general case which 
we leave to future work. We only review the key steps in the derivation of (|17[). since 
the strategy is analogous to the one used in to derive ©• 

We start by recalling the concept of an evaluation representation. The root-of- 
unity representation (|2*9"|) of the finite quantum algebra U q (sl2) can be extended to a 
representation tt^ of the quantum loop algebra U q (sl2) setting 

<(eo) = w^(f), 7r£(/ Q )=™-V(e), <(fc ) = ^(i)~ 2 , 

<(ei) = Tr^(e), = *»(f), ^(h) = ^{tf . (63) 

Here {ej, fa, fei}i=o,i denotes the Chevalley-Serre basis of the quantum loop algebra; see 
e.g. |J| for further details and the conventions used. Employing the coproduct 

A(ei) =ei® l + ki® ei , A(/ 4 ) = £ (g> A:" 1 + 1 ® / i} A(fci) = A* ® A* (64) 

one can build tensor products of representations. In the present context we consider the 
tensor product tt^ ®vr^ of the evaluation representations associated with (|29|) . Without 
loss of generality we can set u = 1. The corresponding representation spaces, which we 
denote by the same symbol, correspond to the auxiliary spaces of the Q-matrices on 
the left hand side of the functional equation (|17j). If the evaluation parameter w is set 
to the special value w = q/ptv the tensor product ir^ <g> k\ becomes decomposable, i.e. 
it contains subrepresentations W\ , W% of the quantum loop algebra giving rise to the 
non-split exact sequence 

Wl ^ 7T{J ® TTi - 1 * W 2 -»• 0, W = q/fJLV . (65) 

Here i : Wi ^> tTw ® n\ is the inclusion and r : 7r^ (g vr^ — ► W2 = vr^ (g ir\/W\ the 
quotient projection. The representations W±,W2 respectively the maps i,t need to be 
determined. This can be achieved by using the intertwiner S(w) : tTw (g 7r^ — > -zr^ (g 7r^ 
detailed in the appendix and exploiting the fact that ker S{q/iiv) = W\. One finds 

Wi = and W 2 = < ® </Wi = tt£, (66) 

with the various parameters given by 

io = q/fJ>v, n = [i-vq-, w = w" = wvq, z = w^iq . (67) 

Here 7T° is the root of unity limit, q 3 — > 1, of the two-dimensional representation of 
U q {sl2) in terms of Pauli matrices and 7r^, the associated evaluation representation of 
the quantum loop algebra. The explicit form of the inclusion and quotient projection 
is given in the appendix. The functional equation (|17|) now follows from the definitions 
(|20|) . (|26|) and the identities 

L» 3 (z/n)L" 23 (zvq 2 )(i®l) = q(z-l)(i®l)L^ q (zfi- 1 i^q)R 23 (zq), 
{ T ®l)L^{z/n)L^{zvq 2 ) = {zq 2 -l)L^{zii- l vq){T ®l) (68) 

which can be verified by explicit calculation using the definitions (|21|) . (|28[l and the 
results (|88|). ()89(l in the appendix. 

Expressing the functional equation (|17j) in terms of the eigenvalues we infer that 
the following ratio 

N' ,,N\t}„{ „N',2N' ,JV\ 



; P B {zq z )P B {zq 
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must be independent of the parameters /x, v. Here we have used the previous results 
(fTB|) and (|5U|) . Consequently, we must have that the ratio 



M^fu/M^ = lim M^=M (70) 

is independent of /u, v. Furthermore, in order that the dependence on /i, v from the 
complete string contribution cancels one is lead to the conclusion (|18|) . Namely, the 
string centre Zj(fi) does either not depend on the parameter [i at all or it just depends 
on it via the simple factor fi 2 . It follows that the ratio simplifies to the expression 

Q,(zq 2 )Q u ( Z qU 2 ) = N z r^ pB{zq 2 )pB{zq) p s{z N'^N = y 

^ifivqyzqv ) 

Inserting this identity into (|17f) and using the previously derived formula (|15|) for the 
eigenvalues of the six-vertex transfer matrix completes the proof of the desired equality 
(HU for N = 3. 

We conclude this section by noting that the result (|18|) now also allows us to derive 
the /Li-dependence of the normalization constant in for real eigenvectors. Employing 
(IKO]) and (J23J setting fi — > fi 1 / 2 , v —* fi one arrives at 

n s 

N»=N^- nB - N ' n ' s , M W * = n^(M 1/2 )/^] JV ' ■ (71) 

Here is simply the number of exact string centres which depend on /i 2 . The above 
identity in particular implies 

M M =M^q-^- nB - N ' n 's (72) 
which fixes the arbitrary sign factor in (|50[) and thus in (|15jl . 



7 Discussion 

In this article we have analysed the eigenvalues of the auxiliary matrices (|26|) at roots 
of unity belonging to the abelian subvariety (JHJ). Let us now summarize what we have 
learnt from their spectra about the eigenvalues and eigenspaces of the six-vertex model 
at roots of 1. 

Our starting point and motivation 1 for employing auxiliary matrices was the ob- 
servation that when q N = 1 the more commonly used approaches such as the coordinate 
space or algebraic Bethe ansatz |15| have serious drawbacks. For example, one de- 
rives from the algebraic Bethe ansatz away from a root of unity the following expression 
for the eigenvalues of the transfer matrix Q2UJI . 

q N ^l: T(z, q )=b(z, q ) M q~^H -* + g»* JJ -* >M (73) 

- - z f (q) z - zf (q) 



.i 



with 



q N ^l: nB = ™-\S*\. (74) 

Here we have set as before z? = e u i q~^. One often finds in the literature that this 
parametrization is used for all real coupling values even though it breaks down when 
q N = 1. This does not mean that the root of unity limit q N — > 1 in (|73j) is ill-defined, 
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it simply requires the explicit knowledge of the Bethe roots. The latter, however, are 
usually not known due to the intricate nature of the Bethe ansatz equations (J2J). What 
can happen in the root of unity limit is that some of the Bethe roots drop out of the 
parametrization (|73j) . As a concrete example consider the spin- zero sector of the M = 6 
chain when q N 1. One finds the three Bethe roots 

M = 6, S z =0: zf = 1, 4 =q-\ z* = q~ 2 

which belong to one of the eigenvalues of the transfer matrix in the four-dimensional 
momentum k = sector. If the limit q 3 — > 1 is taken the products over the Bethe roots 
in (|73[) give simply one and the eigenvalue becomes degenerate with the eigenvalue of 
the pseudo- vacuum consisting of the state where all spins point up (down). 

Since the Bethe roots are not known in general one needs a parametrization of the 
eigenvalues in terms of the finite solutions to the Bethe ansatz equations (J2J) when 
q N = 1. This parametrization we found to be 

q N = i : t(z) = b (z) M q-f+^ n zq2 i f}f + ^- n - n zq ~ 2 ~ zf{q) 



zf(q) ' 1 f = \ z-zf{q) 

where the number of Bethe roots is not fixed by the total spin in contrast to ()74j) . 
Instead, we found the sum rule 



JV 



1 : M — 2noo -2n B = mod N (75) 



for real eigenvectors. Numerical examples for the M = 3, 4, 5, 6, 8 chain and N = 3, 4, 5, 6 
showed so far that it also extends to complex eigenvectors provided / 0. The above 
sum rule also played an important role in the derivation of the Bethe ansatz equations 
(J2J) from the functional equation © . This derivation involved an additional phase factor 
which according to (|75|) is trivial. 

Note that the difference between q N ^ 1 and q N = 1 is not "only" a difference in 
the number of Bethe roots and a change of the phase factors in front of the products in 
(|73|) . In the root of unity limit also the eigenstates of the transfer matrix "re-organize" 
into degenerate eigenspaces across sectors of different spin. The main objective outlined 
in the introduction was to investigate the structure of these degenerate eigenspaces. We 
will now discuss how this information is encoded in the complete strings (|14|) . Recall 
that the complete string contribution in the eigenvalues of the auxiliary matrices (|26|) is 
already fixed by the Bethe root content via the identities ((TH|) and ()19j) . So far we have 
only explained how these results determine the dimension of the degenerate eigenspaces. 
For several examples we will now explicitly see how the complete strings characterize 
the degenerate eigenspaces in terms of irreducible representations of the loop algebra 



7.1 The Drinfeld polynomial and complete strings 

Recall that the loop algebra s/2 has been established as a symmetry of the six-vertex 
model in the commensurate sectors where the total spin is a multiple of the order of the 
root of unity, i.e. 2S Z = OmodiV (S1IZI- in order to make the connection between the 
spectra of the auxiliary matrices and the loop algebra we need first to introduce some 
facts about its representation theory |3U| . 

There are several basis to write down the loop algebra sfa- The most convenient one 
for the present purpose is the mode basis obeying the relations 

h-m+n = [ x m> x n]> [^mj^n] = = ' = 2x m _|_ n , [h m ,h n ] = 0, [l m+1 ,I n ] = [^ ml • (76) 
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The generators {x^,/i m } me z can be successively obtained from the Chevalley-Serre 
basis of the quantum loop algebra Ug es (sl2) at q N = 1 via the quantum Frobenius 
homomorphism Ej (for simplicity we only consider N odd) 

E^^xi, FfUx - EjW^xt, F^^xt t , 2S*/N^ho (77) 
and with the action of the quantum group generators given by [01 Cj^ 

E [ N \q) = F^iq- 1 ) = KE { N \q)ft = K F^ (q' 1 )^ = 

l<mi<---<m JV <Af m^ h m£ h 

Here lEH denotes the spin-reversal operator. As we are only considering spin-chains of 
finite length, all representations of the loop algebra are finite-dimensional and therefore 
highest weight [SO]. That is, there exists a highest weight vector O satisfying 

x^U = 0, h n Q = x^XqQ = x^x~Q, = A n S7, A n G C . (78) 

All finite-dimensional irreducible highest-weight representations are isomorphic to tensor 
products of evaluation representations (HOI- Let tt s : sl2 — ► EndC 2s+1 denote the spin s 
representation of the finite subalgebra sfo = {x^, ho} C sl^- Then define the evaluation 
representation -K s a : sfa — ► EndC 2s+1 by setting 

<(x±)=tt s (x±), K s a {xl x ) = a±V(4), <(/*,)=*'(/*,), aeC. (79) 

The information which evaluation representations are contained in the highest weight 
representation ttq is conveniently encoded in the classical analogue of the Drinfeld poly- 
nomial Pn according to the following correspondence 

n 

7rn^ir£® & Pq(«) = - a jU f s ^ . (80) 

Here all zeroes aj are different. The Drinfeld polynomial can be explicitly calculated 
from the eigenvalues A n of the Cartan elements h n via the following Laurent series 
expansions around u = and u = oo |3U] , 

f> n u" = degPn-u|^, EA_ n u-" = -«^. (81) 

The important observation in connection with the auxiliary matrices constructed in 
is now the following: the Drinfeld polynomials of the highest weight representations 
spanning the degenerate eigenspaces of the transfer matrix coincide with the complete 
string contributions Q19|) appearing in the spectrum of Qu(z) when we identify u = z N . 
That is, we find up to a possible renummeration the identification 

dimvrn = 2 ns and a,- = lim zf(n)~ N . (82) 

At the moment we do not have a general proof of this assertion but we have verified it 
for several examples; see the appendix. We consider one of them in detail to illustrate 
the interplay between the auxiliary matrices and the loop algebra sl2- 

*Here we have used a different convention for the coproduct than in However, one analogously 

proves in this case that the quantum group generators commute with the six-vertex transfer matrix in 
the commensurate sectors. 
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7.1.1 Examples for N = 3 

Consider the spin-chain with M = 6 sites and the primitive root of unity q = exp(2ni / 3) . 
Then the vector with all spins up 

M = 2N = 6: n =t ® T ]= |TT ••• T) (83) 

lies in the commensurate sector S z = OmodiV, where the loop algebra generators are 
defined via (|77|1. The corresponding eigenvalue of the transfer matrix 

T(z,q)\ nn = l + b(z,qf 

is four-fold degenerate with the eigenspace 7tq spanned by 

Here we have indicated via the lower indices the respective spin-sectors. All eigenvec- 
tors have zero momentum. Given the highest weight vector one can now proceed and 
calculate the corresponding Drinfeld polynomial. From the scalar products 



Aq = (o|x^"x n in) = ( ft 



Ai = (n|arfxi|n) = (n 



#(3)^(3) 



^(3)^3) 



^) = I dim7rn = 2, 
= a + + a_ = 20, 

A 2 = (ftlx^xT, \n) = (Q \(x+x^) 2 - \{x^f{xlf\ n) = a 2 + + a 2 _ = 398, ... 
one finds (see also (3T] 1 

P n (u) = (1 - a + u)(l - a-u) with o± = 10 ± 3^ = (10 + SVTl)* 1 . (84) 

Diagonalising the auxiliary matrix Qa(z) in the respective spin sectors one computes 
the following complete string contributions in the subspace of momentum zero, 

S z = +3 : P s (z 3 , /x 6 ) = z 6 - 20zV + ^ = (z 3 - a+/x 6 )(z 3 - a V) 

S z = -3: P s (z 3 ,^) =zV-20zV 3 + l = (z 3 - a+)(z 3 - a_) (85) 

and 

S z = 0: P|(z 3 ,/x 6 ) = z 6 - z 3 (l0(/i 6 + 1) ± ZVli(fj? - 1)) +1 

= (z 3 - a ±f i 6 )(z 3 - a T ) . 

(3) (3) 

Note that {E^'Vl, Ff'fl} are in general not eigenvectors of the auxiliary matrix, but that 
the eigenvectors of Q^{z) are contained in the two-dimensional space spanned by them. 
Taking the limit [i — > 1 and identifying z 3 = u we recover the Drinfeld polynomial (|84j) 
from the complete strings. In this limit the auxiliary matrix becomes degenerate - as 
the representation underlying the definition (|26j) becomes reducible - and E\ F^'n 
are now both eigenvectors of the auxiliary matrix. However, in general we want to keep 
the auxiliary matrix non-degenerate and therefore \x should be chosen different from 
one. 

The above example also nicely confirms the picture outlined in the introduction. 
According to Q18|) there are 2 ns = 4 possible eigenvalues of the auxiliary matrix in the 
degenerate eigenspace of the transfer matrix, all of which we find realized. 

Note that the match between the complete string centres and the evaluation pa- 
rameters is a virtue particular to the auxiliary matrices ()2H|) constructed in pQ. In the 
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context of the six-vertex model the other explicit expression in the literature is Baxter's 
formula (101) in ^U] (which applies to all coupling values 7 6 R but only to the sectors 
of vanishing total spin), 



S Z = : QBa3rter(-8)ai""«M = ex P I 3*7 X) J2 ("m^n ~ "nAJ + \u <*mP' 



Here z = e u q 1 , q = e 1 ^ . Diagonalizing this matrix in the two-dimensional subspace 
{Eq £l,F± Q} of the spin-zero sector we find for each of the two eigenvalues only a 
single complete string with string centres z s = ±1. Thus, for Baxter's auxiliary matrix 
neither the degree of the complete string contribution nor the values of the string centres 
are in agreement with the data obtained from the loop algebra. 

Admittedly, the above example for the M = 6 chain is quite simple and we chose it to 
illustrate the working of the formulas. One might wonder if the identification of complete 
strings and the Drinfeld polynomial also applies when the highest weight vector is a real 
Bethe eigenstate, i.e. when finite Bethe roots are present. We have explicitly worked 
out the following examples with q = exp(2-7ri/3): for the M = 5 chain one finds five 
doublets in the S z = ±3/2 sectors and for the M = 8 chain there are eight quartets in 
the S z = 3, 0, —3 sectors. In all of these cases there is agreement between the complete 
strings (|14|) calculated from the auxiliary matrices (|26|) and the Drinfeld polynomial 
(|8()|) . The results are presented in the appendix. They also show the working of the 
identity Q19[) which yields expressions for the evaluation parameters of the loop algebra 
in terms of Bethe roots. Also the Bethe ansatz equations are recovered by making a 
Laurent series expansion in ()19j) . This shows an intimate link between the Bethe ansatz 
and the representation theory of the loop algebra. 

7.2 Comparison with the eight-vertex model 

We conclude by mentioning that the identity (|19|) coincides with the trigonometric limit 
of a recent conjecture by Fabricius and McCoy [22] on the eigenvalues of Baxter's eight- 
vertex auxiliary matrix constructed in 9 . (Note that this auxiliary matrix is different 
from the one discussed in |1D |I11[ IT2] and p. 3 .) Based on numerical results for the M = 8 
chain and N = 3, 4, 6 they arrive at the elliptic analogue of the identity ()19|) by using 
a functional equation similar in nature to l|17|l. cf. equations (3.12) and (3.10) in |32j . 
However, there are some key differences as Baxter's construction procedure and the one 
used in are not the same. In particular, the dependence on continuous parameters 
analogous to Q respectively (JHJ) is absent in Baxter's matrix j^j. In the context of the 
auxiliary matrices (|26|) we saw that these parameters play a crucial role in breaking 
the loop algebra symmetry and the invariance under spin-reversal. A direct comparison 
between the auxiliary matrices for the eight 1101 1111 I12j and six-vertex model 
is not straightforward: taking the trigonometric limit in the eight-vertex eigenvalues 
requires the knowledge of the explicit dependence of various normalization constants 
and the Bethe roots on the elliptic nome. Further, investigations are needed to clarify 
this point. 
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A The intertwiner for wt <S> k\ with N 



In this section we construct the intertwiner for the following tensor product of evaluation 
representations ir w ig> ix\ with Af = 3. If this intertwiner exists the auxiliary matrices 



Q^(w), Q u (w') must commute, i.e. the conjecture © holds true for N = 3. The defining 
equation of the intertwiner S is given by 

&»«, ® <)A(x) = [« ® <)A op (x)] 5(w), x e C/- ff («J 2 ) . (87) 

Here ir!^ is the evaluation representation (|63[) obtained from ()29|) for A" = 3. The symbols 
A,A op denote the coproduct (jMjl and the opposite coproduct. The latter is obtained 
by permuting the two factors in (|64[). The defining equation ()87|) yields a system of 
algebraic equations for the matrix elements of the intertwiner. As S commutes with 
(nw®TTi)A(ki) it is convenient to decompose the tensor product space into the following 
direct sum 

V = V 1 ®V 2 ®V 3 

where the respective subspaces are spanned by the following basis vectors 

Vi = spanj^o <8 v ,vi <8 v 2 ,v 2 <8 vi}, 
V 2 = sp&n{v (g)Vi,vi(3v ,v 2 (g)v 2 }, 
F 3 = span{i> <8 v 2 , v\ <8 v\, v 2 <8 v } • 

Here V{,i = 0, 1,2 denotes the standard basis in C 3 used in the definition (|29|) of the 
representation tt^. The calculation is cumbersome but straightforward and one finds the 
following solution up to a common normalization factor, 



Hi 



/ 1 






q(wfi—u)(wfj,u—q) 
(wq — fiu)(wq' 2 — fiu) 
w(wfiu—q 2 )(u 2 — q 2 ) 
(wq — /j,u)(wq 2 — fj,u) 



I q 2 (wp,-v) 



5*1 



l 2 



S\ V3 



V 



wq—jxv 


wq—/iv 


w(q—u 2 ) q 


2 (wv—[i) 


wq — fiv 


wq—fiu 








q 2 (w^i—v)(wfi 


-uq 2 ) 


(wq— fiv)(wq 2 


~*%} 


w(w[i—v)(v 2 


-<r) i 


(wq— fjiv)(wq 2 


-aw) 


w 2 (v 2 — q)(v 2 

....\ ...Jl 


-<? 2 ) 





(wfj,v-q)(n 2 -q 2 ) 
(wq — fiu)(wq 2 — fiu) 
q(wv — /j,)(wfii> — q) 
(wq — fiu)(wq 2 — fiu) 




(w^,u—q)(w^iv—q 2 ) 
(wq — fiu)(wq 2 — fiu) / 

(q-^ 2 )(wfj,-u) 
(wq—fiu)(wq 2 — fii') 
[iu(l+q 2 w 2 )+wq(n 2 + l)(u 2 + l) 
(wq— fill) (wq 2 — fj.li) 

w(q— v )(wv— /J,) 
(wq— fiu)(wq 2 — fiu) 



(^-q)(^-q 2 ) 
(wq—fj,u)(wq 2 — fiu) 

(P?-<i )(wv-ji) 
(wq—fiu)(wq 2 — ^iu) 
q 2 (wu — fi)(wu—^tq 2 ) 
(wq—fiv)(wq 2 — fiv) 



We can now use this solution in order to explore the decomposition of the tensor product 
at special values of the evaluation parameter uu. One can explicitly verify that for w = 
qj fj,v the intertwiner has a non-trivial kernel consisting of the following six-dimensional 
space 



ker Si(q//j,u) 
ker S 2 (q/ ^u) 

ker S 3 (q/fjtu) 



span{t>i (8 v 2 , v 2 (8 v±}, 

span{u 2 (8 v 2 , ^—t^-vq (8 «i + vi (8 fo}, 



SP^ Cd-^jf^-l^ O 0^2 + ^2® ^0, V % q _i V Q (g)V 2 +Vi 



This kernel can be identified as a submodule W\ C 7r£ 



tt\ of the quantum loop algebra. 
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A.l The inclusion % : if w , ® tt° z , <^-> W\ C <, ® ix\ 

We define the module W\ simply by stating the inclusion of the basis vectors spanning 
the tensor product 7r^, ® 7r°, into the tensor product -Kw (g) ti\. Denote by {v'j} the basis 

vectors in 7r^, and by {T>|} the basis vector of the two-dimensional representation of 
U q {sl 2 ). Then the inclusion % is defined by linear extension from the following relations 
involving the basis vectors, 



i{y' 2 ® |) = av 2 ®V2, i {v' 2 ® t) = av 2 <S> v±, 
a(Vi<8>T) = 7o ill v i ® v i + V2 <8> v } , i(v[<g) l) = r Y 2 vi®v 2 + ~/ 3 v 2 ®vi, 
i {v'o® t) = Ad u o <8> «i + t>i (8) fo, i(v' ® l) = (3 v ®V2+ i3iVi(g)V\+ (3 2 v 2 ( 3v . 

The coefficients in the above linear combinations are given by 

faV - l)(l - g/iV) 



a 



v([j?q — 1)(^ 2 



2 



a(l - i/ 2 g 2 ) ^(/i 2 - g 2 ) , fl ^ 

7o = M2y2 _ 1 , 7i= g _^ 2 , 72 = (^-1), 

73 = [l a + K<? 2 - "AO ■ 

Acting with the quantum group generators according to (|64[) on the basis vectors in 
the respective tensor products of evaluation representations one verifies that the above 
inclusion is well-defined. 



A. 2 The projection r : vr^ <g> t\\ -> W 2 = vr£, ® ^/Wi 

Having identified the submodule Wi it remains to verify that the quotient space W2 
defines the evaluation representation 71"^,,, as outlined in (|d3|) , (|66j) and (foTj) . This follows 
when identifying the equivalence classes of the following vectors in t\ with the basis 
vectors {v'(} in 7r^„, 

v'o = t(v ®v ), 

v 2 = t(vq (g> V2 — vq 2 v\ ® v\ + v 2 q 2 v 2 (g> vq) . (89) 

This concludes the proof of the decomposition (|65|) , Using the explicit form of the 
inclusion and projection map one is now in the position to proof the functional equation 
<|T7|) as described in the text. 



B Calculation of the Drinfeld polynomial 

We present several examples of calculating the evaluation parameters (|HU|) of the loop 
algebra in the degenerate eigenspaces of the transfer matrix when q = exp(27ri/3). We 
then compare the outcome with the expression (|19|) derived from the complete strings 
of the auxiliary matrices. 



20 



B.l M = 5, S z = ±3/2 

There are in total five doublets for the M = 5 chain in the sectors S z = ±3/2. The 
corresponding highest weight vectors Sl k can be labelled by their momenta and are 
defined as follows, 

5 

n k = £ j nkT (h i) n ITTTTI) , k/K = o, ±2/5, ±4/5 . 

71=1 

Since there are only doublets occurring in this example the corresponding Drinfeld poly- 
nomials Pn fc defined in 1)80)) contain only one factor. For each highest weight vector the 
corresponding evaluation parameter a(k) is calculated using the action of the loop alge- 
bra, 

x+x^n k = [4, + 3qe~ ik + (l + 2q 2 )e- 2ik + (l + 2q)e~ 3lk + 3q 2 e- 4ik \n k 

= a(k)Q k . (90) 

In order to compare this result with the complete strings we may either directly diago- 
nalise the auxiliary matrices (|26|) in the respective spin-sectors or use the identity ()19j) . 
In the latter approach one first solves the Bethe ansatz equation 

and then computes from the Bethe roots the complete strings in the limit — > 1, 

A r (z s_ z s ) _ v {zq2e ~ 1)5 - -A (i i i-io4(^ + g 2 ) J\ 

(92) 

Bethe roots and momenta can be easily matched by taking the limit z — > 1 in 1)15)1 
yielding the following second identity for the evaluation parameter 

a(k) = 10 + Wq 2 /z B - 1/4, e ik = -ql 1 ~ ^ . (93) 

1 — Zb 



Note that the Bethe ansatz equations 1)91)1 are recovered from the Laurent series expan- 
sion in 1)92)1 by setting all coefficients of the terms with powers greater than three equal 
to zero. 

B.2 M = 8, S z = 3,0,-3 

For the M = 8 chain one proceeds similar as in the previous case. One now has eight 
quartets whose highest weight vectors in the S z = 3 sector are again labelled by their 
momenta 

8 

^ = ^e mfc T(i, g niTTTTTTT) . 

n=l 

The degree of Drinfeld polynomial Pn k is now two, i.e. there are two evaluation param- 
eters a± = a±(k) to compute. After some cumbersome computations one obtains 

a+ + a_ = (n k \x+x^\n k ) = 35 + 15ge~ ife + hi^,q 2 e~ 2ik - (5 - 2q)e' 3ik 

+ 6e~ 4ik - (5 - 2q 2 )e- 5lk - 5iV3qe~ m + 15gV 7i/c 
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and 

4a+a_ = (n k \(x+) 2 (xi) 2 \n k ) 

= 4(e~ ik + qe~ 2ik + q 2 e~ Zik + e~ Aik + qe~ 5ik + q 2 e~ 6ik + e" 7ifc ) 2 . 

Again we can compare this result against the complete string by diagonalising the auxil- 
iary matrices or by employing the identity (|19j) . In either case one finds upon matching 
string centres and evaluation parameters the following expression in terms of the Bethe 
roots 

a+ + a_ = 56 + 28q 2 /z B - l/z%, a+a_ = 28 + 56q 2 /z B - 56/4 - 28<? 2 /4 + l/z% . 
Here Bethe roots and momenta are related by the identity 

ik 2 1 ~ z B<f 

e = q . 

1 - ZB 

In order to facilitate the comparison we have summarized the results in the table below. 



momentum 


string centres/evaluation parameters a± 


Bethe root z B jq 2 


k = 


^ (29 ± 3^/93) 

— 


-1 


k = it 


- (83 ± 9>/85) 


1 


k = vr/2 


i (l3(2 + V3)± ^165(7 + 4^3) j 


-2- >/3 


k = -tt/2 


i (l3(2 - V3) ± ^165(7 - 4V3)j 


-2 + v 7 ^ 


k = it/4 


a+ = 38.971..., a_ = 0.0680614... 


1-^-^1(2-2^2) 


k = 371-/4 


a+ = 59.9864..., a_ = 0.615865... 




k = 57r/4 


a+ = 1.62373..., a_ = 0.0166705... 




k = 7vr/4 


a+ = 14.6926..., a_ = 0.0256601... 


1-^ + ^/1(2-2^2) 
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